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Abstract 

Straightforward method for the derivation of linearized version of 
stochastic stability analysis of the nonlinear differential equations is pre- 
sented. Methods for the study of large time behavior of the moments are 
exposed. These general methods are applied to the study of the stochastic 
destabilization of the Langmuir waves in plasma. 



1 INTRODUCTION 



The influence of the external noise on the stability of the linearized dynami- 
cal systems was the subject of numerous studies on finite dimensional models, 
studied both in physical or mathematical literature [T]-|17]. In the following we 
extend part of the formalism to dynamical systems described by partial differ- 
ential or integro-differential stochastic equations. 

The main motivation of this study is related to the sensitivity of the onset 
to unstable regime, under the effect of random noise, observed in hydrodynamic 
turbulence pQ, the effect of stochastic destabilization of the stable regime of 
simplest Hamiltonian systems, like linear undamped LC system with random 
capacitance [2], [3]. The main dynamical system of interest for the controlled 
thermonuclear fusion studies with magnetic confinement are described the ki- 
netic equations. Before the study of these more complete models, we expose 
the main physical and mathematical principles on a more simple model: of 
Langmuir waves in turbulent background. 
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The large part of the mathematical aspects exposed here, at least those 
applied to finite dimensional systems, can be found in the mathematical litera- 
ture. An example of application to large discrete systems, perturbed by space 
and temporal white noise, can be found in the works [10j-|13j. 

The derivation of the equations for the stochastic linearized system from the 
stochastic nonlinear one as well as the application to Langmuir waves, is new. 

1.1 Physical model 

The main class of problems studied in this article are related to the influence of 
the external noise on the linear stability of steady state of general, non linear 
physical systems. We suppose that the system is described by a set of general 
nonlinear partial differential, or integro-differential equations, that contains a 
stochastic perturbation, modelled by a random field. 

In all of this part of the study the random field will be modelled by temporal 
white noise with general spatial correlations. 

Instead of establishing a general abstract framework, in this article we illus- 
trate the method by an example. The first interesting physical model is non the 
linear hyperbolic evolution equation, that describe the propagation of nonlinear 
waves in plasma, under the effects of the external noise n u (x, t). 

= ~ n^(x,t),n u (*,t)) (1) 

Here F(u, v) is a nonlinear function, having near u = v — the linearized 
form F(u, v) = au + /3v + O(u 2 +v 2 ),withu>0. We are interested in the 

u,v— >0 

linear stability of the solution in the deterministic case $(x, i) = 0, under the 
effect of the noise n(x, t) , represented in fact by a random field whose statistical 
properties are supposed to be known. In particular we suppose that the random 
field is homogenous and stationary. 

We illustrate now the general recipe to reduce the linear stochastic stability 
analysis of non linear systems to standard form. Similar to the deterministic 
case we expand F($,n) in the variable $ up to order 2 term to obtain 

F($,n u ) =a$ + Af w (x,O* + a w (x,t)+0($ 2 ) (2) 

where a > is constant, and the multiplicative and additive noise term are 
given exactly by 



M w (x,t) 



d[F($,n u (x,t)) - F{$,0)] 
<9$ 



(3) 



a u (x,t) = F(0,n u (x,t)) (4) 

Neglecting the 0($ 2 ) term and introducing the variable P w (x, t) = 
we obtain the linear system of stochastic differential equation of the first order 
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form 



d*„(x,i) = P^dt 

dP u (x,t) = (A + a)$dt + dW u (-x.,t)$> + a u (x,t)dt 
dW w (x,t) = M u (x,t)dt 



(5) 
(6) 
(7) 



Without loss of generality we can consider the case 



In the following we will approximate M u (x, t) by a temporal white noise, con- 
sequently we have 



Here C(x, x') is the spatial correlation function of the random field, related 
to the random multiplicative terms. By using L p estimates like in the work 
[S], it can be proven a result similar to [BJ, that under general conditions the 
additive term does not provide exponential destabilization of the moments of 
the random fields $ w , P u . Consequently we consider in the continuation only 
the homogenous system, while we are interested only in the large time behavior 
of the linearized system Eqs.([5][7|). 

The previous system of stochastic differential equation can be formally rewrit- 
ten as general linear stochastic differential equation in a suitable Banach space 
as follows 



Similar linearization procedure can be performed for the system of drift 
kinetic equations, describing various species of charged particles (q a , m a ) , 1 < 
a < ^4in a strong constant magnetic field B. The starting nonlinear system of 
equations for the functions / Q (x, v,t) can be written as follows 



E w \W u (x,t)] =0 
E w [W w (x,i)W r u (x',t / )]=n 



min(i,OC(x,x') 



dY u (t)= Adt + dB u (t) Y u (t) 



(8) 



where 




dfg(x,v,t) 
dt 




dfa q a „ df a 



E x B 



B 2 



a 
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The linearization of the equations starts with the ansatz 

/ (x,v,t) = 5f u>a + [n (x) + <5n w (x,t)]exp [-v 2 (j8 Q (x) + %,(x,i))] 

around the equilibrium density and inverse temperature profiles no(x), /3q(x) 
. The stochastic effects are induced by the temperature and density fluctuations 
(5/3 w (x, t), <5n w (x, t) . By performing the same linearization procedure we obtain 
a system of linear stochastic differential equations of the form Eq.© where 

/*/ u , 1 (x,v,t)\ 

v _ 5/u,,2(x,V,*) 

\<y/w,^(x,V,t)/ 



1.2 Mathematical background 

We will study now the correct formulation of the stability problem of the solu- 
tions of Eq.©, and describe some general methods to solve these problems. 

Consider a general linear evolution equation in a complex Banach space V 
with the norm || .| . 



dt 



= Ax(t); x(f) e V; (9) 



where A is a linear operator with dense domain in V and range in V.We 
will denote by V* the dual space of V, where the transpose A T of A acts. 

An equivalent form, written in components ( or reducing our study to the 
discretized form of the partial differential of integro-differential equations) is 

d ^=A]^{t) (10) 

where summation on repeated indices is understood. 
Typical examples of interest will be the linearized equations used in the 
deterministic linear stability studies. 

In the framework of the deterministic linear stability studies, in the finite di- 
mensional clear answer to the stability problem is given by the eigenvalue 
of A with largest real part. In general case we denote by sp(A) the spectrum of 
A and ^ 

A = max|Re(z)|z € sp(A)\ (11) 

The real number A is the largest Liapunov exponent. Then for all e > we 
have for any solution of Eqs.(ji?l fTUj) 

hm - 1 = (12) 

t->ooexp[(A + £)i] 

If A is known, then the question about the stability of null solution in Eqs.Q 
fTUj) is solved. 
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Consider now the stochastic generalization of Eqs.® [TU|), by adding a ran- 
dom, time dependent a linear term, for the sake of simplicity modelled by an 
operator -valued Brownian motion . In this end, we consider (CI, A, FtjP) a 
filtered probability space and denote by B u (t) , with co € fi, an operator- valued 
stochastic process adapted to the filtration J- t . In the our case B u (t) is an oper- 
ator valued Brownian motion, or equivalently _B w (i)x is an V valued Brownian 
motion [IB] for all x G V. 

Remark 1 We will use everywhere the ltd formalism. 
We have the following stochastic evolution equation, 

dx u (t)= Adt + dB u {t) x^t) (13) 

or in components ( we use the general relativity index conventions) 

dxt(t) = [A)dt + plJw^t)] xi(t) (14) 
1 < i,j < n; 1 < a < A (15) 



where a are constants and m° (t) are independent standard Brownian mo- 
tions 

E w [< (t)v& (*')] = min(i, t>)6 a < b (16) 



We have the matrix- valued Brownian motionsi? w (t), with components B* u (t), 



having the correlation tensor Cj' n 

^ [B}Jt)B™Jt>)} = mxn(t,t')C}'™ 

A 

/~ii,m \ ^ i m 

j,n / j rj,arn,a 



(17) 
(18) 



o=l 



To this correlation tensor we associate a linear operator C acting in the 
tensor product space V <g> V, so an equivalent form of Eqs. (fTTI IT5|) is 



E^ 



{ [B u (i)x] <g> [^(tOyl } = min(t, t')C (x«y);x«yeV«V (19) 



An equivalent form of the Eq. pl?]) that will be useful is 



{ dB u (i)x <g> dB u {t)y }=C(x<g)y)cft 



(20) 



Remark 2 Observe immediately that exactly due to the ltd definition, if we 
denote Yi(t) = Ex w (f), i/ien we obtain from Eg. U3\) 



dY^t) 
dt 



AYi(t) 



(21) 



so apparently the effect of the noise is averaged out. For higher order moments 
the equations will contain also the correlation effects of the noise. 
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Nevertheless, we can try generalize Eq. (fT2| as follows: 

Problem 3 Find the minimal value of X p such that for p > 1 and all e > we 
have 

Mg#ig = (22) 
t^oc exp [(Ap + e) i] 

Simplest soluble one-dimensional cases, as well as more general results [I], 
suggest that in this case the exponent A is dependent on the exponent p in a 
non-trivial way: in the case of systems that in deterministic case are linearly 
stable, it is possible that we have Eq . (|2"2")) for small noise intensity and small 
values of p we have A p < 0, hence stability, but surely for sufficiently large p we 
have Ap > . 

Example 4 Consider the one dimensional model 

dx u {t) = (-adt + pdw u {t))x u (t) 

in the ltd formalism, where Wu,(t) is the standard Wiener process (E u w u (t) — 
0; E w [w tl j(t)] 2 =t), a, a are constants. The general solution is 

x u (t) = x(0) exp[-(a + p 2 /2)t + aw u (t)] 

. Results 

[E w M*)H 1/P = K0)|exp(Apt) 
with A p = —a + (p — l)p 2 /2 . 

Conclusion 5 In the stochastic case the linear stability is determined by a scale, 
Ap of Liapunov exponents 

The function X p was effectively computed only in some special cases ([B], 
0). 

We will see that excepting the small noise, perturbative calculations, in 
the physically interesting cases ( partial differential equations, Vlasov-Poisson 
systems) computation of A p for even values of p is possibly only at the expense 
of solving eigenvalue problems with spatial dimension increased by a factor p . 

In this moment there are no general method to compute A p , for all real 
values of p, even for the case when A or B u (t) are 2x2 matrices. So we 
reformulate the Problem( [3]) in a more weak sense as follows 

Problem 6 For a given degree p, find X p such that for all e > and all mono- 
mial m p (x) of degree p > l(i.e. m p (x) = x^x 1 " 2 ■ ■ ■ x ln with i\ + • • • + i n = p ) . 
the case p = \ is trivial, see Remark^) ). In the components Xi we have 

lim E - kg (X " M = (23) 
t^ocexp[p(X p + e)t] 

In this end we will obtain closed sets of deterministic equations for the 
moments 

E u) [x%(t)xZ(t)---xt(t)] (24) 
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2 The deterministic equations for the moments. 
2.1 The evolution of the mean values 

For simplicity consider first the finite dimensional case, denote in general x = 
(x 1 , x 2 , ...x n ).We have the following 

Remark 7 (Backward Kolmogorov equation) Consider a general stochastic pro- 
cess x tL ,(t), described by the following SDE and initial conditions 

dxl(t) = V l (^{t))dt + < (x w (t))d<(t); 1 < i < n; 1 < a < A (25) 
^ [ < (*)] = 0; E w [ < (t)wt (t')] = 6 a > b min (t,f) (26) 
x Wiy (0)=x (27) 

For for any differentiable function F(x) we have (for the sake of simplicity of 
the notations, the subscript uj will be omitted) 



BF 

dF(^(t)) = — [V\^{t))dt + al{^{t))dwl{t)] 



1 d 2 F 



D^ixjdt 



2 dx^dx 1 
A 

D^(x) = ^4(x)^ (x) 

a=l 



We denote 

Then we have 

Mpit) 
dt 



M F (t) =E U [F^it))] 



dF{^{t)) 1 #F . 

V J + 7;7r-7r-D (*)) 



dx l 



2 dx^dx 1 



(28) 

(29) 
(30) 
(31) 



It follows from Eqs. ([3"Tl |2"9")) that in the case when the SDE (|2"5j) is linear, 
it is possible to obtain a closed set of linear differential equations when F(x) 
belongs to a finite dimensional subspace of homogenous polynomials. 

We consider now the case when Eg. ([231) is a linear SDE, like Eq. (fT4")) . So 



A)x> 



(32) 
(33) 



We obtain from Eqs.(|S51[551 [THJ 

D^(x) = C^ n x m x n (34) 
We obtain from Eqs. (|3Tl l32l f34|) the basic equation for further development 



M F (t) 
dt 



9F(xu(t)) f , . . I d 2 F ■ ■ 

Q xl A i X uW + 2 dxj Q x i L rn,n *u W X u (t) 



(35) 
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2.2 Linear equations for the moments of order N. 
2.2.1 CaseiV=l. 

Consider, for more clarity reasons, in Eq.® F(x ) = J x k , where k is fixed 
integer index and denote in general (t) = [ x£(f)] , as in Remark @. 

Recall, x k is a contravariant vector component, not power. 

Then we obtain Eq. (j2~Tj) . or 

jyfe 

1 =4 y i(*) ( 36 ) 



dt 
dYi(t) 



AYi(i) (37) 



The corresponding stability problem is the classical, deterministic 

Problem 8 Compute the eigenvalues, left and right eigenvectors in the eigen- 
value problems 

Au = Am; u € V; 

Av = A T v : v 6 V*; v is left eigenvector 
2.2.2 Case N = 2 

Consider in Eq.JHJ F(x ) d = x k x l and denote F 2 fc ''(f) = E w [ a£ (t)a£, (*)] . Then 
we obtain 

dY k ' 1 

= A k Y^{t) + A]Y^(t) + C k m l n Y™- n {t) (38) 

In condensed notation, taking into account that to Y^' l {t) we can associate 
a tensor Y 2 (i) £ V ® V and to C^|„ the operator C (see Eq. (f2"U)l ) we have the 
equivalent form of Eq. 



dY 2 (t) 



A®1 + 1&A+C 



Y a (t) (39) 



2.2.3 Case AT = 3 

Similarly we consider _F(x ) = x k x l x i and denote Kf^f) = E W [ x k (t)xl{t)xl(t)] 
. Then we obtain 

= A kyj,u +A )Y*' j ' i +A)Y*> l < j + (40) 

< ~'m,n / 3 '^m,!!^ T"°m,n- r 3 

In order to write in a more compact form, we introduce Y 3 (i)eV®V®V 
with components Y^' 1 ' 1 and the operators Ci j2 , Ci^, and C 23 , acting in 
V <g> V (g) V according to the last terms in Eq. (j4T)|) , and we obtain 
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= (A®l®i + l(g)A(g)l+l®l(g)A+ (41) 

Ci,2 + Ci,3 + C2,3)Y3(£) 



2.2.4 The general m case 

m 

With the choice F(x ) d = ~[[x k i and notation Y**'"' ' km (t) = E u 

3=1 

from Eq. ([35|) results 



3 = 1 



jyfci,-- ,fc m 

r-^- = Ayi l fe l -,U +i4 fc ay fej,-,*« + AyhrJ+ (42) 

at J J J 

Esik a ..kb-w-ki_,---k a _ 1 ,j,— ,fe m 
'in 

Similarly we have the compact form 



[A<g)(lJ + I<8u4<8 f lj + (43) 

9(m-l) _ . „ 

<E>A + }_^C a . b }Y m (t) 

a<b 

de J- U m (C)Y m (t) (44) 

where the action of the operator C a ,b m the tensor product space V® m is 
defined according to the corresponding term in Eq. (j42| . Then the large time 
behavior of Y m (i) is dominated by the term exp(A m i)Z m where A m , in the 
finite dimensional case, is the eigenvalue with largest real part in 

AZ m = U m (C)Z m (45) 

For application of the perturbative analysis in the limit of weak noise inten- 
sity (in Eq.(j43| the contribution of C a ,b is small) it is useful to remark that, 
Eq. (j45"l) can be put in the form 



AY r , 



(M m +<5M m )Y m (46) 



where the unperturbed term is 

i(m— l) ^ ^ ®(m— 2) 



M m = A® (l J + f 1J + (47) 

• •• + (!] «>4 
and the perturbation is 

SM m = ]Tc M (48) 

a<b 



9 



Remark 9 If the deterministic stability problem\$\is solved then the spectrum 
(including eigenvalues), , left and right eigenvector of M m from Eq. ^jTtycan 
be computed immediately. 



2.2.5 Perturbative approach. 

Perturbation of general, not necessary self adjoint operators. These 
results are generalizations of the perturbations methods from quantum mechan- 
ics to the case when the operators are not necessary self adjoint. 

The formalism exposed above is useful to derive perturbative methods, to 
solve the eigenvalue problems like from Eg . (|45[) .Consider a general complex 
linear vector space U Suppose that the deterministic linear stability analysis 
was performed, so we know the eigenvalues, v kl as well as the right eigenvectors 
u,R,fc € U and left eigenvectors v^.fc e U* of the operator M . The left 
eigenvectors, in matrix notations, can be associated to the eigenvectors of the 
transposed M T . Denote the linear function on U associated to x* £ U* by 
x — > (x* , x) . The transposed M T is defined by 

(x* , MxJ = ( A? T x* , (49) 

The spectrum of M and M T identical in very general case. In particular 
in finite dimensions it results from the property of determinants. Much of our 
results are explained in the finite dimensional case and we do not discuss the 
specific problems related to continuos spectrum and unbounded operators . 



Simple eigenvalue So to each eigenvalue Vk corresponds the right eigen- 
vectors ur,*; £ U, and the left eigenvectors VL,k G U* . We have 

M\i R , k = v k U-R,k (50) 
M v Ljk = is k v L ,k (51) 

Proposition 10 The eigenvalues v% + 5vk of the perturbed operator M + SM 
can be computed in first approximation order by 

(vL,k,SMu Rt k) 
Sv k = ±- (52) 

(v L ,fc,U fl , fc ) 

Proof. The perturbed right eigenvector will be represented as ur.^+Jur^. Results 
(M + SM^ (u R . k + Su Rtk ) = {v k + 5v k ) (u R . k + Su Rtk ) (53) 
We use now [501 and drop the second order terms. We obtain 

M6u R ,k + SMu Rnk = v k Su Rjk + 6i/ k u Rt k (54) 
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By left multiplication with u^^ results 

v L , k ,M5u R}k ) + (v£ ifc ,<SMu fl ,fc) = v k (v Ljk ,5u Rjk ) + 5v k (v Ljk , u Rtk ) (55) 



By using Eqs. (l49l [STj) , the first term in l.h.s. in Eq. ([55| can be rewrit- 
ten as (v Lik , M5u Rlk ^ = (^M T VL,k,5u R>k ) = v k (vL,k,8u Rtk ) .Combined with 
Eq. ([551) we obtain Eq. (l52"1) . that completes the proof. ■ 

Multiple eigenvalue, semisimple (without Jordan blocks) Must of 
physical cases from this category are related to some discrete or continuos sym- 
metry group, and the explicit symmetry breaking due to asymmetric perturba- 
tion ( Like Zeeman or Stark effects). Let v k an eigenvalue of multiplicity m. 
The subspace of left and right eigenvectors has the same dimension. Let VL,k,a 
respectively Ul k,a > with 1 < a < m the set of left and right eigenvectors 

Mu Rjk , a = VkU R ,k,a (56) 

M T VL, fe , a = v k v L<k<a (57) 
Define the m x m matrices F and G by matrix elements as follows 

F a ,b = (vL,k,a>Mu R>k ,b) (58) 

G a ,b = (vL,fe,a,U_R !fc!& ) (59) 

Denote by bv k a = z k with 1 < a < m , the sequence of m roots z k of the 
secular equation in z 



det 

We have the following 



F-zG 



= (60) 



Proposition 11 The multiplicity, m semisimple eigenvalue v k under the effect 
of perturbation splits into m eigenvalues 

v k -> v k + Su k>a (61) 

Proof. Combination of the previous proof and from the degenerate case from 
quantum mechanical textbooks. ■ 

2.2.6 The case of diagonal noise 

Our equations will be simplified in the particular case, when the noise in SDE (fT4l) 
has the form 

P\ a = fyi (62) 

In this case 

C*>™ = 6]6™\q> m (63) 
^2' m = X>WC (64) 
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The Eq.flU]) became 

dYb'-' K 



A ki Y j,k 2 ,-,k m +A k 2Y ki,j,...,k m +A }»yferj + ( 65 ) 



fa J m 3 m J 

a<b 



3 Langmuir waves in turbulent background 
3.1 Statement of the problem 

The evolution of the electrostatic potential in Langmuir waves in homogenous 
finite temperature plasma, after rescaling x — > \JiVthermaV<- is given by the Klein- 
Gordon equation 



9^(x,t) 



= A<p - mT<i> (66) 



dt 2 

2 

2 2 n e e 



e m e 



In the case of perturbation of the background electron density n e — > n e 
const C( x j^) by a temporal white noise with 

E[C(x,i)] =0 
E[C(x,t)C(x',0]=£(t-i')C(x,x') 



formally can be written as 
9 2 0(x,i) 



= A x 0(x,t)- [m 2 + C(x,t)] 



dt 2 

Or in a more rigorous form, by introducing the Brownian motion io(x, t) 
oh that formally £(x, t) = dw ^^ j we obtain the system of first order PSDE 

# w (x,i) = p u (x,t)dt (67) 
dp w (x, t) = -K^udt + 4>dw u) {y., t) (68) 

where -ftT x = to 2 — A x and E w [u> w (x, t)w u (x', t')] = min(i, t')C(x, x') . 
Denote 

F^(xi,x 2 ,i) = E w [$w(xi,i)$a>(x 2 ,*)] 
5^(xi,x 2 ,i) = E„ [p w (x 1 ,t)$ w (x 2 ,t)] 
y^ 1 p(xi,x 2) *) = E^ [$ UJ (xi,i)p w (x2,t)] 

y p) p(xi,X 2 ,f) = E w \p u ,{x. 1 ,t)p u (-K2,t)] 
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By the same procedure we obtain the following system of linear equation for 
the second order correlation functions 

rfy^(xi,x 2 ,f) = y ^ (xi ^ ^ t) + y ^ (xi ^ ^ ? t) (g9) 
dr p ,^(xi,x 2 ,t) = _^ r ^ (xi ^ i} + y ^ (xi (7Q) 

d^, p (xi,x 2 ,t) = rpp(xijX2ji) _ ^ 2 y^( Xl)X2)t) (7i) 

dy p , p (xi,x 2 ,t) = _^y 0j) ( Xl)X2)t) _^y pi0(xi)X2)t) (72) 

+ C(xi, x 2 ) y0,0(xi,x 2 ,f) 

In order to study the t — > oo behavior of the solutions of Eqs. (|69H72p . we 
consider the particular class of solutions 

y a ,5(xi,x 2 ,t) = exp(At)Z„ j6 (xi,x 2 );a, 6 € {p, </>} (73) 

and after simple algebra we find the following problem: 

Problem 12 Find A with maximal real part from the generalized eigenvalue 
problem 

A 4 + 2A 2 (2m 2 - A X1 - A X2 ) + (A Xl - A X2 ) 2 - 2AC (x 1; x 2 )] Z^(x 1; x 2 ) = 

(74) 

The boundary conditions are obtained in a straightforward way from the 
boundary conditions of the original deterministic problem. 

3.1.1 Alternative notations 

A more transparent way is given by a compact notation. We will denote 

f*A />(x,t)\ 

^(x,*y 

Define the operator valued 2x2 matrices L x and J as follows: 

L* = | I T 

-K x 

0^ 

1 6, 

or by components 
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Then Eas. (|671 155]l became 

d*(x, t) = L^dt + JVdw (75) 

We use the notation 

F il>i2 (xi,x 2 ,t) =E w [tf illU (xi,t)¥ il>ttf (xa,t)]; l<*i,*a<2 (76) 
Respectively Eqs. (|69M72j) became 



ri^i,t 2 ( x i> x 2;£) _ Xj 2 

^ ^ ii ,ji r n M + L i 2 , 32 r h ,32 {") 

+ J nJiJi2,]2 Y n,32 C ( x i> x 2) ; 1 < H,l2, h, h < 2 

Consider the stability problem of the moments Yi lj i 2 (xi, x 2 , t) . From the 
representation 

Fi,i 2 ( x i> x 2,£) =exp(iujt)Z ilii2 (x 1 ,x 2 ) (78) 
and Eq. (|77[) the following eigenvalue problem results 

iuZ iui3 (xi,x 2 ) = L^ h Z jlti2 + Lfl h Z ilih (79) 

+ J! </i 2 , j 2 ^il J2 C ( X l I X 2 ) 

In analogy to the Eq. (|43[) similar equations for the higher order correlation 
functions can be obtained. 

3.2 Particular solutions of the eigenvalue equation (I74L 
791). 



3.2.1 Zero noise limit 

This is useful for small noise limit perturbative calculations. 
In line with the classical stability analysis we will denote 



e k = y/m 2 + k 2 (80) 
A = iu (81) 

The classical dispersion relation, for the Langmuir waves is 

w k = ±£k (82) 

When C — we have the unperturbed eigenvectors in Eq. (f79f of the form 

Zi 1 ,i a (xi,x 2 ) = v iu i 2 exp [i (kixi + k 2 x 2 )] (83) 
v n,i2 = Ui t u i2 (84) 

By using Eqs. (f79l l83l l84| after simple algebra we find that the generaliza- 
tion of the classical dispersion relation, for the case of two point correlation 
function, from Eq. ([82)l is, 

w kl ,k 2 = ±£ki ± £k 2 (85) 
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3.3 

3.3.1 Complete spatial correlations ( |17| ^ . 

This part is included only for the purpose of clarification the application of 
formalism. Can be considered as a limiting case , of the results exposed in 
subsection (j3.3.2j) 

Consider the simplest soluble case C (xj., X2) = er 2 = const . In this case 
the ansatz from Eq. (|83|) can be used and from Eq.([79|) results (see Appendix 

U Eq.CCED) 

telA* + 2iuj<j2 = ( £ k 1 + 4 2 - ^ 2 ) 2 

or by using Eq. ([8lj) 

A 4 + 2A 2 (4 i+ 4 2 ) + (4 i -4 2 ) 2 = 2Aa 2 (86) 

From Eq.( l86|) results that if |ki| — |ki2 1 is sufficiently small then we have 
two real solutions with A > , so the random multiplicative perturbation of 
the background always produce parametric destabilization of the 2 point corre- 
lation functions Y < p_ < p(xi,x 2 ,t) = E w [$ Cl} (xi,t)$ u (x2,t)]. Due to the fact that 
Y^ ^(xi,X2,£) is symmetric with respect permutations Xi±^X2 , the modes with 
to the modes (now labelled by a couple (k^k^)) with |ki| ~ |k 2 | has non- 
vanishing contributions. For |kj| = |k 2 | we have for the dominating mode 

A 3 + 4A4 = 2 a 2 

For high |k| the the exponential growth is dominated by A x a 2 / (2E 2 ,) + 

|k|— >oo 

0(l/s^) . The most dangerous modes are those with for |ki| s» |ka| ~ 0. For 
small a we have A max X cr 2 /(2m 2 ) + C(cr 6 ) 



3.3.2 Homogenous system 

Suppose that the random field is homogenous, so C(xi, x 2 ) = C (xi — x 2 ) , 
with C {— x) = C(x). We introduce the "center of mass" and relative distance 
coordinates (xi + X2) /2,r = xi — X2 and use the ansatz in Eq. ([7H) 

Z , (xi, x 2 ) = exp [ik (x x + x 2 ) /2] *(r) (87) 

The resulting generalized eigenvalue problem is 

A 4 + 2A 2 (2m 2 + k 2 /2 - 2A X ) + 4 (ikV x ) 2 - 2AC (x)l *(r) = 



Temporal and spatial white noise. |17j we study now the opposite 
case studied previously in subsubsection (13.3.11) . Consider the case one di- 
mensional version of the previous model, with Dirac-Delta correlation function 
C(xi, x 2 ) = <7 2 <5(xi — x 2 ). It is known that the bound state problem in the 
Schrodinger equation in higher dimension with delta function potential give 
mathematically inconsistent results 22 . 
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The one dimensional version of Eq. (1551) is 



A(X, k)i/>(x) - B(X, k)i>"(x) - 2\<t 2 6(x)^(x) = (89) 

A(X,k) = A 4 +4A 2 m 2 + X 2 k 2 (90) 
B(X,k) = 4(A 2 + k 2 ) (91) 

From Eq.(??) results the continuity and jump conditions for the restriction 
of tp±(x) to the domains x > respectively x < 0. 

]imip-(x) = lim^+fx) (92) 



B(X,k) 



lim — lim^_(a;) 



x/0 



+ 2Acr 2 V>+(0) = 



(93) 



The class of physically admissible (because A is complex) solutions are of the 
form 



ip+{x) = exp(— ax) 
tj)—{x) = cxp(ax) 
Re(a) > 



By straightforward calculations from Eqs 

B(X,k)a 2 = A(X,k) 
2B(X,k)a = 2Xa 2 

a 4 = 4(A 2 + fc 2 )(A 2 +4m 2 



results 



And we obtain 



k 2 ) 



(94) 
(95) 

(96) 
(97) 

(98) 



From Eqs. ([93 ??) results that for any root of Eq.® with Re(A) > , the 
condition from Eq. ([95|) is fulfilled. 

It is clear that now there is a threshold: when cr 4 < 4fc 2 (4m 2 + k 2 ) then A 
is imaginary. For cr 4 > 4k 2 (Am 2 + k 2 ) there is a real positive A, so the system 
is unstable. We remark an important mechanism: for low noise intensity first 
the large wavelength structures are destabilized. So by this mechanism there is 
transfer of energy from short to the long wavelength. The most sensitive to the 
destabilization are the mode with small k. 



Positive correlations . Consider the solutions of the Eq. 
C(x) > and lim C(x) = 0. 



in the case 



Case k = We are interested in the case A ^ 0, so the generalized eigen- 
value problem can be reformulated as follows. 



A x (r) 



1 

2A 



C(x) 



A 2 

<Jr(r) = ~( m 2 + _ )$( r ) (99) 
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When < A < +00 , , because C > 0, the lowest eigenvalue E\ of the 
eigenvalue problem- A x *(r) + [ -^C(x)] *o,A(r) = Ex^a,\{v) 



- A x (r) 



*o,A(r) = ^ A *o,A(r) 



-(rn 2 + ^) = E x 



(100) 
(101) 



is a monotone increasing continuos function, from —00 to 0. 

Indeed, the term —j^C (x) is like an attractive potential well in a ground 
state problem in the Eq. (|100p . formally identical to Schrodinger equation. For 
A — > it is very deep potential well, can be proven rigorously that E\ lim E\ = 

—00 . When A increases then there is a critical value \ cr it > such that 
the bound state problem from Eq. (jl00j) has no more solutions . So we have 
lim E\ = 0. (Previous asymptotic estimates can be verified by approxi- 

mating the potential well with a rectangular one, at least in the case of fast 
correlation decay) . 

Results that there is a one and only one value of A > such that E\ = 
— (to 2 + 4-) (because the r.h.s. is decreasing, from — m 2 to —00) . 

Conclusion 13 It follows that there exists at least an eigenvalue A > in the 
generalized eigenvalue problem from Eg. i99\) , so at least the mode k = is 
destabilized for arbitrary small noise intensity. 



Case |k| — > 00. In general, Eq. (|go|) . by suitable choice of the coordinate 
system (k ||Oz), can be reformulated as follows. 




*o,A(r) 



#o,A(r) = ^A*o,A(r) 

(102) 
(103) 



,2 k 2 A 2 
~(- 2 + T + T ) = £a 



where Aj_ 



D.v- 



a 2 

9y 2 



Because of the new term (-|) , the ground 
state energy has a lower bound. 

Be denoting with 5 the typical width of the maximum of C (x) near x = 0, 
respectively c = C(0), for large k and small A we have the estimate E\ > 
0{— S ) . So for fixed large values of k there is no real A that satisfies 
Ea. p5l l . 

Conclusion 14 . For k > there is a threshold for the intensity of the multi- 
plicative noise such that the mode remains stable below this threshold. 
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4 Conclusions. 

Similar to the simplest one dimensional case exposed in the works [2], [3], the 
Langmuir waves in a turbulent background are destabilized. The energy of the 
fluctuations related to the multiplicative noise is transferred with an exponential 
rate to the long wavelength fluctuations preferentially. In contradistinction with 
one dimensional case, there is a threshold in the noise intensity for each mode, 
such that below this threshold the mode remains stable. 

5 Appendix 

Denote the tensor with components Vi 1: i 2 from Eq. (j53"|) by v . We use the 
notation Eq.(??) and define 2x2 matrices i?(k),a as follows 



and by 1 the 2x2 unit matrix. The eigenvalue equation (|79[) with ansatz 
Eq.([531) became 




(104) 



(105) 



(106) 



(107) 



icov =Bv 



(108) 



where 



B = £(ki) <g) T +1 ®E(k 2 ) + cr 2 aig>a 



(109) 



From Eq. (llOSp results 



o; 2 v =i3 2 V 



(110) 



We will use the following matrix identities, that can be verified easily 








(111) 
(112) 
(113) 



Ui + u 2 



1 



From Eqs.(Hni [ml HI3 EES]) results 




(114) 
(115) 
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and with the notation 



y 



the eigenvalue equation (|110|) became 

Cv = yv 

with 

C = 2S(ki) <g> E(k 2 ) +(T 2 (l <g> a + a (8 1 



It follows 

pil 2 

G v = y v 
and from Eq. (|118[ 1111111131 ) we obtain 



C 2 - 44 x 4 2 + 2<7 2 



S(ki) ® 1 + 1 ®E(k 2 ) + cr 2 a®a 



or compared with Eq. (|109[) 



.2 2 



So, from Eqs. (|1191 112ip we obtain 



4e 2 i£ 2 2 + 2a 2 B 



v = 2/ v 



or, combined with Eqs. (|108l II we obtain the relation 
44 i e 2 2 +2 lW a 2 = ((4 i+ 4 2 )-^ 2 ) 2 



(116) 

(117) 
(118) 
(119) 

(120) 
(121) 
(122) 
(123) 
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